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Abstract 

Four subcritical response (SR) methods were 
evaluated fOr on-line use in transonic wind- 
tunnel tests where the flutter model is excited 
solely by airstream turbulence. The methods 
were: Randomdec, Power-Spectral -Density, Peak- 

Hold, and Cross-Spectrum. Subcritical response 
data were obtained during tests in the Langley 
Transonic Dynamics Tunnel of a cantilevered 
flutter model wing. The test procedure was to 
maintain a constant Mach number (M) and increase 
the dynamic pressure (q) in incremental steps. 

The test Mach numbers were 0.65, 0.75, 0.82, 

0.90, and 1.15. The four methods provided damp- 
ing trends by which the flutter mode could be 
tracked and extrapolated to a flutter-onset q. 

A "hard" flutter point was obtained at M = 0.82. 
The Peak-Hold and Cross-Spectrum methods gave 
reliable results and could be most readily used 
for on-line testing. At M = 0.82, a p-k analy- 
sis predicted the same flutter mode as the exper- 
iment but a 6-percent lower flutter q. At the 
subcritical dynamic pressures, calculated damping 
values were appreciably lower than measured data. 

Nomenclature 

A amplitude of spectrum peak (Fig. 2) 

f,f_ frequency of vibration mode corresponding 
to peak in spectrum, Hz (Fig. 2) 

g damping coefficient (Fig. 2) 

M Mach number 

N number of wave cycles (Fig. 2) 

q dynamic pressure 

^CALC flutter q calculated at M = 0.82 

qp flutter q predicted by damping trend 

qp^ average qp 

q,p q of last (highest-q) point used in 
damping trend extrapolation 

qpi flutter q measured at M = 0.82 

SR subcritical response 

Xq X|. amplitude of Randomdec signature at 
’ . start and after N cycles, respectively 

(Fig. 2) 

a(t), h(t) time history of model torsion and 
bending motion response, respectively, 
(Fig. 2) 

Af frequency bandwidth at half-power level 
of PSD peak, Hz (Fig. 2) 


I. Introduction 

Flutter tests in wind tunnels, are made to 
define the flutter characteristics of models of 
complete aircraft or their components. Such 
testing necessarily entails the risk of model 
damage due to the destructive oscillations that 
can occur at flutter. This risk is reduced if 
the flutter onset can be predicted from the 
model behavior in the subcritical region below 
the flutter boundary. To accomplish this, the 
vibration modes critical to flutter must be 
identified and the damping and frequencies of 
these modes measured and tracked as the test 
conditions are varied until a damping trend can 
be reliably extrapolated to a flutter condition 
of zero damping. Reference 1 describes several 
methods used in wind-tunnel and airplane flight 
flutter tests for extracting these modal quan- 
tities from subcritical structural responses. 

The present paper reports the results of an 
evaluation of four subcritical response (SR) 
methods for on-line use in transonic wind- 
tunnel flutter tests where the model response 
is due solely to airstream turbulence. 

The four SR methods evaluated were Peak- 
Hold, Power-Spectral -Density (PSD), Randomdec, 
and Cross-Spectrum. The first three of these 
have been used in previous wind-tunnel tests. 2 
The Cross-Spectrum method is new and attempts 
to measure quantitatively the coupling between 
the bending and torsion motions that occur as 
flutter is approached. The response data used 
in the evaluation were acquired during a 
typical flutter test conducted in the Langley 
Transonic Dynamics Tunnel (TDT), The purpose 
of this test was to determine the transonic 
flutter characteristics of a research model 
wing. The evaluation was made to establish the 
SR method(s) most applicable during flutter 
tests of this type and was based primarily on 
the practical use and capability of each method 
to provide subcritical damping trends from 
which the flutter-onset boundary could bo 
accurately estimated. 

The flutter research model used in this 
study represented an aeroelastically tailored 
wing proposed for a fighter airplane. This model 
was designed and built by General Dynamics/Fort 
Worth under Air Force contract and tested in 
cooperation with the National Aeronautics and 
Space Administration. The aeroelastic tailoring 
aspects of the study are not presented in th'is 
report. The semi span wing model was tested 
cantilever-mounted as shown in figure 1, Model 
subcritical response data were obtained during 
flutter tests at M = 0,65, 0,75, 0.82, 0.90, and 
1.15. A "hard" flutter point was obtained pur- 
posely at M = 0.82, 

Post-test analyses of the magnetic tape 
records of the test model -response data were 
made using the four SR methods. Procedures 
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were Implemented so that the data were reduced 
and Interpreted within time limits consistent 
with on-line test use. For comparison with 
experimental data, the model flutter charac- 
teristics and flutter-mode damping trend were 
calculated at M = 0.82 using a "p-k" 
analysis. 3 Presented herein are the results 
of thU study. 

II. Description of Methods 

An overview of the data reduction procedure 
for each SR method is presented in Figure 2. 

All of these methods could be implemented on 
the data acquisition system (DAS) of the TOT 
which includes a Xerox Sigma 5 digital 
computer. 4 The TDT DAS was used only with 
the Randomdec method. High-speed frequency 
analyzers were used with the other three SR 
methods. All four methods can process the time 
histories of the model responses on-line, i.e., 
while the test is in progress, with the present 
equipment. An arbitrary time limit of about 
1.5 minutes was set to acquire, reduce, and 
interpret the response data; this permitted a 
data sampling time of from 30 to 60 seconds. 

In near-flutter regions, adequate damping mea- 
surements can be obtained in much less time. 

The time history of the torsional motion 
transducer response a(t) was processed in all 
SR methods (Fig. 2). The Cross-Spectrum method 
employed the bending motion transducer response 
h(t) also. The torsional response was selec- 
ted as the primary response signal because it 
gave the best damping trends. 

Peak-Hold Method 

In the Peak-Hold method, the response time 
history is processed and converted to a fre- 
quency spectrum by a Spectral Dynamics SD330A 
Spectroscope-Real Time Analyzer. This analyzer 
was operated in the Peak-Hold mode with the 
spectrum frequency range set from 0 to 100 Hz. 
For this frequency range, the response time 
history is processed in time samples of 2.5 
seconds duration with some overlapping of the 
samples. Displayed on the analyzer scope was a 
continuously updated frequency spectrum of the 
response amplitude. Specifically, the value in 
each of the 250 frequency windows of the spec- 
trum was the maximum value of the root-mean- 
square response amplitude measured in any one 
time sample during the time the analyzer was 
operated (typically, 60 seconds). 

In the resulting spectrum the frequency of 
each response peak corresponded to that of a 
vibration mode. It has been found, e. g. , in 
refs. 5 and 2, that tracking the inverse of the 
response amplitude value for a given mode gives 
a good indication of the damping trend to a 
flutter condition (see Fig. 2). 

Cross-Spectrum Method 

The Cross -Spectrum method provides a quan- 
titative measure of the coupling between the 
torsion and bending motions that occur as flut- 
ter is approached. This method is an attempt to 
quantify what had been regarded in many past 
tests as an indication of flutter, namely, the 


regularity and duration of the Lissajous pat- 
tern formed from measurements of the wing bend- 
ing and torsional motions. Time histories of 
the bending and torsional motions of. the model 
were input separately to a Hewlett-Packard 
Digital Signal Analyzer, Model HP5420A. The 
analyzer was set for a cross -spectrum measure- 
ment in the frequency range from 9 to 59 Hz. 

In this setting, the responses are processed in 
ensembles of time history samples of 5 seconds 
duration with some overlapping of the 
ensembles. Displayed on the analyzer scope was 
a continuously updated cross-spectrum which was 
obtained by stable ensemble averaging. A 
50-second time history of response data was 
processed during which 60 ensembles were 
averaged to form the final cross-spectrum. 

In a cross-spectrum, each response peak 
could usually be associated by frequency with a 
vibration mode. As was done in the Peak-Hold 
method, it was assumed that the amplitude of a 
response peak was Indicative of the damping in 
that mode, i.e., the greater the response, the 
lower the damping. A logarithmic relationship 
between the damping and the cross-spectrum 
amplitude was found to give generally consis- 
tent damping trends and was therefore used in 
this study (see Fig. 2). 

PSD Method 

The PSD bandwidth method was implemented in 
a manner similar to the Cross-Spectrum method 
except that the HP5420A analyzer was set for an 
auto-spectrum measurement. A value of the 
damping coefficient g was calculated from the 
response peak bandwidth at the half -power point 
as shown in Figure 2. 

Randomdec Method 

For the Randomdec method^, a 30-second 
time history of the response signal was digi- 
tized by the TDT data acquisition system. A 
fast Fourier transform (FFT) of the time his- 
tory was performed to acquire a frequency spec- 
trum of the signal. From the frequency spec- 
trum, a peak response was chosen for further 
analysis of the time history data. The time 
history was passed through a recursive digital 
band-pass filter with a center frequency 
corresponding to the frequency of the chosen 
peak response. The filtered data were then 
ensemble-averaged to form a Randomdec signa- 
ture. Damping values of the signature time 
history were calculated by using a linear 
least-squares fit of the logarithmic values of 
the signature peaks (see Fig. 2). 

III. Evaluation Procedure 
Evaluation Criteria 

Each SR method was evaluated for its capa- 
bility during on-line testing to provide the 
following: 

1. Accurate flutter-onset prediction. 

2. Consistent damping trend. 

3. Early flutter-onset prediction as 
flutter approached. 

4. Continuous on-line damping indication. 



5. Actual damping value. 

6. Ease and simplicity of use. 

Model and Wind-Tunnel Test 

The model (Fig. 1) was a 0.25-size, dynam- 
ically and elastically scaled version of a 
fighter airplane wing which was geometrically 
similar to an F-16 airplane wing. This wing 
was aeroelastically tailored for washin aero- 
dynamic characteristics. Some model physical 
properties, including measured vibration char- 
acteristics, are presented in Fig. 3. Strain 
gages were mounted near the wing root and 
oriented to measure the bending and torsional 
motions of the model. 

The model was tested in Freon* in the 
Langley Transonic Dynamics Tunnel. The model 
wing was cantilever-mounted to a turntable on 
the wind-tunnel sidewall. An aerodynamic 
fairing, resembling a fuselage half-body, 
allowed the wing root to be located outside the 
wall boundary layer (Fig. 1). The model was 
tested at angles of attack that gave a near 
minimum aerodynamic load. An attempt was made 
to define a flutter boundary at Mach numbers of 
0.65, 0.75, 0.82, 0.90, and 1.15. In these 
tests, M was held constant while q was 
slowly Increased in Increments until a flutter 
condition was considered reasonably predicted 
using the Peak-Hold and PSD methods. The model 
was excited solely by the airstream turbulence. 
Model strain-gage response data were acquired 
at dwell points of constant q. At M = 0.82, 
the model was purposely tested to a "hard" 
flutter point. 

Flutter Analysis 

A flutter analysis using the p-k method3 
was conducted for M = 0.82 to provide 
subcritical damping data for comparison with 
the test results. In the analysis, the 
generalized masses were determined using 
calculated mass distributions and the measured 
frequencies and mode shapes of the first five 
vibration modes. These mass distributions were 
obtained from a finite-element mathematical 
simulation of the model wing structure and had 
been checked and corrected (if necessary) to 
agree with mass measurements of the model 
component parts and the complete model. The 
off-diagonal, cross-coupling terms in the 
generalized masses were found to be relatively 
small (thus indicating the measured mode shapes 
were reasonably orthogonal) and were neglected 
in the flutter analysis. 

The unsteady aerodynamic terms in the 
analysis were generated by subsonic lifting- 
surface (kernel function) theory. 7 For the 
analysis, thirty-six aerodynamic collocation 
points were distributed over the wing surface, 
with six chordwise points located at each of 
six spanwise stations. The analysis provided 
data at matched wind-tunnel test conditions of 
velocity and density. 


*Freon: Registered trademark of E. I. 
duPont de Nemours & Co., Inc. 


IV. Presentation of Results 

The results of the flutter analysis are presen- 
ted in Figure 4. The analysis predicted the 
model would flutter at M = 0.82 in the funda- 
mental bending mode (mode 1) at a dynamic pres- 
sure of 8.47 kPa (177 psf) and a frequency of 20 
Hz. In the test at M = 0.82, a "hard" flutter 
point was measured for the model at a dynamic 
pressure of 9.00 kPa (188 psf) and a frequency 
of 20.8 Hz (Table I). 

Model response data obtained at the test dwell 
points were processed using the data reduction 
procedures previously described for each SR 
method. With each SR method, a damping indicator 
value and associated frequency were determined 
for the 20-Hz flutter mode. Damping indicator 
values and frequencies were also determined for 
other identifiable vibration modes, but these 
data are not presented because they were not the 
primary indicators of flutter. Plots of the 
damping values and associated frequencies against 
dynamic pressure are shown for each test M in 
Figures 5 through 9. A discussion of these 
figures follow. 

A practical problem in using these SR methods 
is how best to define a damping trend for extrap- 
olation to flutter onset. In this study, all 
data points that exceed a rough scatter band 
about an obvious trend were neglected. Also 
neglected were damping values at the lower dyna- 
mic pressures that did not significantly affect 
the flutter-prediction trend. The neglected 
damping values are identified by the flagged 
symbols in Figures 5 through 9. To define a 
damping trend, an attempt was made to fair 
through the data points a second-order curve 
using a least-squares fit. However, if an unde- 
sirable curve shape resulted such as a concave-up 
curvature that gave an unconservative flutter 
dynamic pressure estimate, it was replaced by a 
linear least-squares fit. Because this curve 
fitting requires some judgment and considerable 
time, it is suggested that on-line test damping 
trends and flutter-onset extrapolations be 
checked by thorough post -test data reduction. 

The dynamic pressures for flutter onset (qp) 

obtained from extrapolating each damping trend to 
a zero damping value are compiled in Table I and 
plotted in Figure 10(a). For comparison the 
calculated and measured flutter dynamic pressures 
at M = 0.82 are indicated in Figure 10(a). A 
measure of the scatter in qp values predicted 

by the different methods from the average qp^ 

is included in Figure 10(b) in t erms of an 
equivalent flutter speed ratio /q 

The damping (g) values calculated by the p-k 
analysis at M = 0.82 are compared with those 
measured by the PSD and Randomdec methods in 
Figure 11. For the comparison, the g values 
are plotted against dynamic pressures that have 
been normalized by either calculated or measured 
flutter dynamic pressure of Rp|. 

To determine the lowest q level at which 
flutter onset can be reliably predicted, the 
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damping data at M = 0.82 were reworked as 
follows. From each set of damping data points, 
the damping value for the high6st-q point was 
dropped and the curve fitting and extrapolation 
redone to obtain a new qp value. This 

represented the qp that would be predicted 

if the test had been stopped at the last, 
highest-q point (q|_p) of the remaining damping 

data. This procedure was repeated until there 
were only two damping points left for 
extrapolation. The resulting qp values are 

shown in Figure 12. In a similar manner, 
flutter-onset predictions using calculated 
damping values (mode 1 in Fig. 4) were determined 
and are included in Figure 12. Lines of equal 
q values = <lLp) ore drawn on these 

plots. Normalized plots of these data are 
presented in the upper portion of Figure 12. The 
dashed, cross-hatched lines on the plots define a 
variation of +5 percent in q. 

V. Discussion of Results 
Analysis-Test Comparison 

At M = 0.82, the calculated and measured 
flutter point data are in good agreement. The 
p-k analysis predicted a 6-percent lower flut- 
ter q and a 4-percent lower flutter frequency 
than did the experiment. This correlation con- 
firmed that the flutter mode had been correctly 
identified and tracked in the tests. However, at 
dynamic pressures near the flutter point (Fig. 
11), the calculated g values are appreciably 
lower than those measured either by the PSD or 
Randomdec methods. The reason for this discrep- 
ancy is not known. These data were verified by 
thorough checks of both the analysis and experi- 
mental procedures. 

The calculated damping trend with q is 
approximately linear in the proximity of the 
flutter point (Figs. 4 and 11). The damping 
trend for the Randomdec method appears linear 
(Fig. 7), and the trend for the PSD method has 
only a slight curvature. The damping trends 
for the Peak-Hold and Cross -Spectrum methods 
have pronounced curvatures. Thus, although the 
g levels of the PSD and Randomdec methods do 
not agree well with calculated values, they do 
exhibit roughly similar damping trend shapes. 

Accurate Flutter-Onset Prediction 

The average (^p^) value is considered to 

represent the actual flutter-onset q value at 
each test Mach number except 0.82 and 0.9. At 
M = 0.82, the test "hard" flutter point is 
available for comparison. At M = 0.9, the model 
was buffeting considerably and there was concern 
that the measured damping trends were not truly 
indicative of flutter proximity; therefore, the 
test was ended at a rather low q for fear of 
losing the model before the other Mach numbers 
were explored. Unfortunately, the wind-tunnel 
schedule did not allow further testing at this 
M. Consequently, at M = 0.9 the SR results 
have by far the greatest scatter and the actual 
flutter-onset q is least certain. The 


experimental flutter boundary drawn in Figure 
10(a) is the best estimate of the, boundary 
shape with the dashed lines indicating the 
uncertain portions. 

All four SR methods predicted nearly the 
same qp at M = 0.65 and M = 0.82 

(Fig. 10). At the other Mach numbers, there 
were varying amounts of scatter in the qp 

values predicted by the different methods. As 
might be expected, the closer the model was 
tested to flutter, the closer grouped were the 
predictions. At M = 0,82, comparison with the 
"hard" flutter q shows that the qp 

values predicted by the SR methods (Table I) 
range from about 1 percent low to 2.7 percent 
high (PSD method), with the average predicted 
value (qp^) nearly exactly equal to the 

measured test point (q|^). At the other 

Mach numbers, no method appears to be conclu- 
sively superior in predicting accurately the 
qp^ value (Fig. 10(b)). 

It was found that a qp value was 

dependent on the number of data points used to 
establish the damping trend. For example, at 
M = 0.82, where damping data were obtained very 
close to the flutter point, the qp values 

were slightly different for damping trends that 
used either all the data points or only the 
last five highest q points (see Table I). In 
the present extrapolations, all the data points 
that formed a consistent damping trend were 
used because it was felt that they would give 
more accurate results. 

Consistent Damping Trend 

At all test Mach numbers, the damping 
trends for all four SR methods were fairly con- 
sistent over the test dynamic pressure range 
(Figs. 5-9) with the PSD and Randomdec having 
the most scatter in the damping values. The 
damping trends for the Peak-Hold and Cross- 
Spectrum methods can usually be identified and 
tracked at lower dynamic pressures. This 
could save testing time by allowing the test 
conditions to be changed more rapidly than when 
there were no indications of damping levels. 

It is of interest to note that no SR method had 
the same damping-trend shape at all test Mach 
numbers. Nor did all SR methods have the same 
damping-trend shape at any one Mach number. 

Early Flutter Prediction 

All four SR methods begin to predict the 
flutter-onset q with reasonably accuracy at 
roughly the same test q level. This is 
demonstrated in Figure 12 which shows that to 
predict the flutter onset at M = 0.82 within 
+5 percent, the test must be conducted to a 
dynamic pressure at least within about 7 to 10 
percent of the actual flutter point q^. 

Although the Peak-Hold and Cross -Spectrum 
methods provide fairly consistent damping 
trends at the lower dynamic pressures (Fig. 7), 
the qp values predicted from these low-q 
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data are not particularly accurate (Fig. 12). 
Table I shows that at all Mach numbers except 
0.90, the tests were conducted to q levels 
that were within at least 9 percent of the 
value. Therefore, these qp^^ 

values are believed to be reasonably good 
estimates. 

The qp vs, qj^p plots of Figure 

12 provide indications of when to terminate a 
test at a particular M, It can be seen that 
as flutter is approached, the qp values 

tend to converge and approach an intersection 
with the Pp = qj^p line at the actual 

value. Thus, a test could be ended 

when the Pp values have converged 

to an acceptably constant value and/or the 
intersection with the qp = q^^p line 

is reasonably well defined. Both the Peak-Hold 
and Cross -Spectrum methods appear to do this at 
lower q levels than the other SR methods. It 
is suggested that, in using these SR methods, 
similar plots be made as a routine part of the 
data reduction. 

Continuous Damping Indication 

As presently implemented, the Peak-Hold 
and Cross-Spectrum methods can provide nearly 
continuous damping trend data more readily than 
the other two methods. This is because the PSD 
and Randomdec methods required more operator 
involvement in the data reduction process. It 
is very helpful to have a measure of damping 
while changing conditions during a test. 

Actual Damping Value 

Only the PSD and Randomdec methods provide 
actual g values. However, the g values ob- 
tained by these two methods do not agree with 
each other at most of the Mach numbers (Figs. 
5-9) nor, as mentioned before, do they agree 
with the g values calculated by analyses at 
M = 0.82. The reason for this discrepancy is 
not known. In defense of the experimental 
methods, they do provide consistent damping 
trends that give good qp values. The 

Randomdec method can extract damping values at 
low-q test conditions where the PSD method is 
inapplicable. The accuracy of a PSD measured 
damping value becomes poorer as the damping 
level increases because the response PSD peak 
shape from which the damping value is deter- 
mined becomes irregular and less well defined. 
The Randomdec signature is usually well defined 
over at least several wave cycles (Fig. 2); 
therefore, there is more confidence in the mea- 
sured damping value. At the low-q, high damp- 
ing test points, the accuracy of the damping 
measurements for both of these SR methods would 
probably be improved if longer than the present 
30 to 60 seconds of response time history was 
processed. 

Ease/Simplicity of Use 

The Peak-Hold and Cross -Spectrum methods 
were the simplest and easiest to use. The 


Randomdec and PSD methods require some judg- 
ments and occasional data reworking, 

VI, Concluding Remarks. 

Four subcritical response (SR) methods, 
namely, the Peak-Hold, Cross -Spectrum, Power- 
Spectral -Density (PSD), and Randomdec methods, 
were evaluated for on-line prediction of flut- 
ter onset in a wind-tunnel test where the model 
was excited solely by the airstream turbu- 
lence. The present evaluation was limited to 
the application of these SR methods in a single 
but typical flutter test of a cantilevered 
flutter model wing. The test Mach numbers were 
0,65, 0.75, 0,82, 0.90, and 1.15. More experi- 
ence using these methods with different model 
configurations and test conditions is needed. 

All four SR methods provided flutter-mode 
damping trends from which the flutter-onset 
dynamic pressure could be reliably predicted. 
The Peak-Hold and Cross -Spectrum methods are 
recommended for on-line test use because they 
provide damping trends more easily and rapidly 
than the other SR methods. 

For M = 0.82, a p-k flutter analysis was 
made using lifting surface (kernel -function) 
theory. Comparisons of the analytical results 
were made with the experimental data for the 
"hard" flutter point and SR damping values. 

The analysis predicted the same flutter mode as 
the experiment but a 6-percent lower flutter 
dynamic pressure. At the subcritical dynamic 
pressures near flutter, the calculated damping 
values were appreciably lower than the experil 
mental data. 

It was found that to obtain a reliable 
flutter-onset prediction from SR damping 
trends, the test must be conducted to within 7 
to 10 percent of the actual flutter dynamic 
pressure. It is suggested that on-line damping 
trends and flutter-onset predictions be checked 
by thorough post -test data reduction, 
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Table I.- Summary of test results. 





qp predicted by 





Mach 

Number 

^LP 

Peak-Hold 

Cross- 

PSD 

Randomdec 

‘Ifa 


‘’m 

number 

of 


Spectrum 




Data 

Points 

kPa 

kPa 

kPa 

kPa 

kPa 

kPa 

‘’fa 

kPa 

0.65 

all 

10. 43^ 

10.68 

10.49 

10.52 

10.39 

10.52 

0.99 

- 

.75 

all 

9.00^ 

9.24 

10.44 

9.20 

9.62 

9.63 

.93 


.82 

all 

8.86 

8.96 

8.97 

9.24 

8.92 

9.03 

.98 

9.00^ 


last 5^^ 

8.86 

9.10 

8.95 

9.10 

8.95 

9.03 

.98 


.90 

all 

6.80^ 

9.15 

10.05 

7.34 

7.09 

8.41 

.81 

- 

1.15 

all 

11. 92^ 

13.45 

13.45 

12.60 

12.69 

13.04 

.91 

- 


♦Maximum test q. 

♦♦Extrapolated only five highest-q points. 
^"Hard" flutter frequency was 20.8 Hz. 
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Fig. 1 Model mounted in wind tunnel. 


Fig. 2 Implementation of suticritical response 
methods. 
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(a) MODEL DIMENSIONS IN INCHES. 
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Fig. 3 Model physical properties. 
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Fig. 6 Damping trends measured by four SR methods. 
M = 0.75. 
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Fig. 4 Damping trends predicted by analysis. 
M = 0.82. 
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Fig. 7 Damping trends measured by four SR methods. ■ 
M = 0.82. 
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Fig. 8 Damping trends measured by four SR methods. 
M = 0.90. 
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Fig. 5 Damping trends measured by four SR methods, 
M = 0.65. 



q. kPa q, kPa q, kPa q, kPa 

la) PEAK-HOLD lb) CROSS-SPECTRUM (c) PSD Id) RANDOMDEC 

Fig. 9 Damping trends measured by four SR methods. 
M = 1.15. 
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(a) Flutter boundary predicted by SR methods, 
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(b) Variation in flutter speed predictions. 

Fig. 10 Predicted flutter boundaries and 
comparisons. 
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Fig. 11 Comparison of measured and calculated 
damping values. M = 0.82. 
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Fig. 12 Flutter onset predictions using different 
damping-trend end points. M = 0,82. 
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